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ABSTRACT 

We present a general method for testing gamma-ray burst (GRB) jet struc- 
ture and carry out a comprehensive analysis about the prevalent jet structure 
models. According to the jet angular energy distribution, we can not only derive 
the expected distribution of the GRB isotropic-equivalent energy release for any 
possible jet structure, but also obtain a two-dimensional distribution including 
redshift z. By using the Kolmogorov-Smirnov test we compare the predicted 
distribution with the observed sample, and find that the power-law structured 
jet model is most consistent with the current sample and that the uniform jet 
model is also plausible. However, this conclusion is tentative because of the small 
size and the inhomogeneity of this sample. Future observations (e.g., Swift) will 
provide a larger and less biased sample for us to make a robust conclusion by 
using the procedure proposed in this paper. 

Subject headings: gamma rays: bursts — gamma rays: observations — ISM: jets 
and outflows — methods: statistical 



1. INTRODUCTION 

Growing evidence such as the achromatic breaks in afterglow light curves of ~ 20 
gamma-ray bursts (GRBs) and the observation of polarized emission (Covino et al. 1999; 
Wijers et al. 1999) suggests that GRBs are produced by collimated jets. Numerous models 
for the jet configuration have been proposed due to its importance in understanding the 
burst progenitor and GRB event rates. One leading model is the uniform jet one (e.g., 
Rhoads 1997; Sari, Piran, & Halpern 1999; Frail et al. 2001; Dai & Cheng 2001 ), in which 
the angular energy distribution within the jet is uniform while the nearly constant total 
energy is collimated into different opening angles. Another is the structured jet model, in 
which the jet energy and structure are approximately identical for all GRBs, but the energy 
per solid angle e{0) varies as a function of angle from the jet axis within the structured jet 
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(Meszaros, Rees, & Wijers 1998; Dai & Gou 2001; Rossi, Lazzati, & Rees 2002; Zhang & 
Meszaros 2002). Furthermore, although the specific distribution of jet opening angles in the 
uniform jet model is uncertain, several kinds of distribution such as the power-law, Gaussian 
and exponential function have been suggested. In this paper we discuss these possible models 
specifically. 

Since a few leading models of jet structure can all explain the observed achromatic 
breaks in afterglow light curves, some authors attempted to test the geometrical configu- 
ration by using other observations. For example, the jet structure could be constrained 
by considering the luminosity function (e.g., Firmani et al. 2004; Lloyd- Ronning, Dai, & 
Zhang 2004; Guetta, Piran, & Waxman 2005). The logiV — logP plot was also used to 
constrain the jet opening angle distribution for the uniform jet and the star forming rate 
(Lin et al. 2004; Guetta, Piran, & Waxman 2005). Perna, Sari, & Frail (2003) predicted the 
observed one-dimensional (ID) distribution of viewing angles, n{6) = dn/d9, in the power- 
law structured jet model. They found that the predicted distribution could fit the current 
sample of 16 bursts with modelled angles and known redshifts. Considering the same model, 
Nakar, Granot, & Guetta (2004) carried out an analysis about the two-dimensional (2D) 
distribution n(z,6) = d 2 n/dzd6. However, they found that the 2D prediction shows a very 
poor agreement with the same sample, and the hypothesis that the data are drawn from 
this model is rejected at the 99% confidence level by a 2D Kolmogorov-Smirnov (K-S) test. 
Based on the standard energy reservoir of GRB jets and the relationship between the cosmic 
rest frame GRB spectral peak energy and the isotropic gamma-ray energy, i.e., -E pca k oc E^, 
Liang, Wu, & Dai (2004) simulated a sample including 10 6 GRBs to test the GRB proba- 
bility distribution as a function of viewing angle and redshift, and supported the power- law 
structured jet model. Dai & Zhang (2005) summarized the above tests and performed a 
global test on the quasi-universal Gaussian-like structured jet. Together with their previous 
tests with the observed jet break angle vs. isotropic energy (E iso — 9j) and observed peak 
energy vs. fluence relations (E°^l k vs. fluence), Zhang et al. (2004) concluded that the 
current GRB data are generally consistent with the Gaussian-like jet model. These tests 
can only tentatively support certain possible configuration of the GRB jet, partly because 
the size of the sample use in the test is too small (but see Cui, Liang, & Lu 2005, who 
recently adopted the BATSE sample). In this paper, we propose a new method to test the 
GRB jet structure by using the sample of bursts with measured fluences and redshifts. By 
considering three jet structure models, i.e., the power- law structured jet model, the Gaussian 
structured jet model, and the uniform jet model, we work out the theoretical distribution 
of burst isotropic-equivalent energy release in both ID distribution N(E) = dN/dE and 
2D distribution N(E,z) = d 2 N/dEdz (hereafter, we use E instead of E iso for brevity), and 
compare them with the observed distribution of E derived from the sample with measured 
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fluences and redshifts. 

This paper is organized as follows: In §2.1 we present several prevalent models for 
jet structure and the probability distribution of the isotropic-equivalent energy under these 
models, while in §2.2 we describe our methodology in detail. In §3 we discuss what each 
model predicts for the observed distribution (in ID and 2D) of GRB isotropic-equivalent 
energy release and assess how our numerical results are compared with the current sample 
by using K-S test. We summarize our results and draw our main conclusions in §4. 

2. THEORY 

2.1. Models for jet structure 

Since the energy release in an explosion plays an important role in determining the 
progenitors of GRBs and the physics of their central engine, some authors have carried out 
plenty of analyses about the gamma-ray energy release distribution. Within the uniform jet 
model, Frail et al. (2001) found that the gamma-ray energy release after beaming correction 
is narrowly clustered around 5 x 10 50 ergs by fitting a sample of 16 well observed GRBs. 
Using a larger sample, Bloom, Frail, & Kulkarni (2003) obtained a similar result that the 
gamma-ray energies are well clustered around a median value of 1.33 x 10 51 ergs. The gamma- 
ray energy release should become larger if the gamma-ray radiation efficiency is considered 
(Xu & Dai 2004). Panaitescu & Kumar (2001, 2002) modelled the broadband emission of 
several GRB afterglows and obtained a similar value of the kinetic energy. Furthermore, 
Wu, Dai, & Liang (2004) derived a jet break time-flux density relationship and constrained 
some physical parameters of gamma-ray burst afterglows. The jet break relation supports 
the "standard candle" hypothesis of the afterglows by an entirely different approach. We 
therefore adopt a standard value for GRB jet energy Ej = 10 51 ergs throughout this paper. 
According to the presently measured isotropic-equivalent energy, we adopt the upper limit 
of the isotropic-equivalent energy as -E max = 10 55 ergs. 

According to the angular energy distribution for any plausible jet model, one can obtain 
the isotropic-equivalent energy as a function of 9. The isotropic-equivalent energy, E, is 
defined as E(9) = Aire(9) ) where 9 is the angle from the jet axis and e(9) is the energy per 
solid angle. 

In the power-law structured jet (PLSJ) model, all GRBs are supposed to have a similar 
jet profile, and their different observed properties are due to different observer's viewing 
angles. Under the power-law jet scenario, the angular distribution of the isotropic-equivalent 
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energy is 

where Ie determines the sharpness of the transition from the core to wing, and represents 
the profile of the jet wing. The opening angle of the core C is introduced to avoid the 
divergence of the total energy. E c is defined as E c = E(0), and equals to the maximum 
isotropic-equivalent energy -E max - Below we assume Ie = +00 for simplicity. If all bursts 
were observable, then in a structured jet model, the probability for a viewing angle is 
P(0)dO OdO. From the relation P(0)dO = P(E)dE, we obtain the normalized probability 
density at a given isotropic-equivalent energy E, 

o Q2 / rp \ l+2/fc E 

F(£) ^Sfe(i) • ^W k ^<E<E c ). (2) 

In the Gaussian structured jet (GSJ) model, the isotropic-equivalent energy distribution 
that varies with the angle from the jet axis is (Zhang & Meszaros 2002) 

_ e 2 

E(0) = E c e *%, (0 < < tt/2), (3) 

where E c = E(0), and C is the characteristic jet angular width. According to the same 
derivation as in the power-law structured jet model, we have the normalized probability at 
a given isotropic-equivalent energy E, 

P(E) = ^ 6 -l, (E c e- 2 ^ 2 < E < E c ). (4) 

7T Z E 

Similarly we assume E c = 10 55 ergs. 

In the uniform jet (UJ) model (e.g., Lamb, Donaghy, & Graziani 2005), all GRBs 
produce jets with different jet opening angles 0j. The structure of a uniform jet reads 

m ~U 3 <0<n/2. (5) 

Assuming f(E) is the probability that a GRB releases the isotropic energy between E and E+ 
dE, the probability that we observe a specified E is given by P(E) = f(E)(l — cos 0j). Using 
the Frail's relation, Ej = E{1 — cos 0j) ~ constant, we obtain the probability distribution of 
the isotropic energy 

P(E)nf(E)^, (Ej < E < E max ), (6) 
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Next we turn to discuss the specified form of f(E). As for the GRB opening angle 
distribution, Lin et al. (2004) studied three cases in which the GRB opening angle 0j follows 
a power-law, exponential, or Gaussian distribution. Here we also consider these three cases. 
We assume that the first probability distribution of isotropic-equivalent energy follows a 
power-law function with index (, 



f(E)= pl+ ^ C pl+c iS C , (EjKEKE^). (7) 

-CTnax III a 



It is worth noting that f(E) should be normalized, i.e. /J' max f(E)dE = 1. The second 
probability distribution for f(E) is an exponential function with parameter A, 

\p XE 3 

f{E) = T —^r^ e"-, ( Ej <E<E max ), (8) 

which has also been normalized. Finally, if f(E) is a Gaussian distribution with the mean 
energy E and standard scatter a, we have 

f(E) = -^T^, (Ej < E < E max ) : (9) 
where a is the normalization factor satisfying J i ^ max f(E)dE = 1. 



2.2. Detections 

For a power-law photon spectrum, the photon luminosity in the triggering band v% — u u 
is given by 

E a - 2 (vi\~ a 1 - (vu/vi) 1 '" 



LphM Thv { a - 1 U J 1 - (^ 2 /^i) 2 - Q ' (10) 
where a is the photon spectrum index, and v\ and v u are the lower and upper limiting 
frequencies of a GRB detector respectively. For the BATSE triggering channel, v\ = 50 keV 
and v u = 300 keV. The E is the isotropic-equivalent 7-ray energy in a range of v\ — 0.1 keV 
to U2 = 10 4 keV at the GRB's cosmological rest frame. The mean spectral peak energy for 
GRBs is (E p ) ~ 250 keV (Preece et al. 2000). According to the spectral fit, we take a ~ 1 
for the frequency range v x < v < v u in the following calculations (Band et al. 1993). Then 
the photon luminosity can be simplified as 

W*j = Th{ £_ - } ln (^) . (ii) 
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where T is an "effective" duration that the burst would have if its energy output is constant 
at the peak value rather than the highly variable one. 

The detector can be triggered if the photon flux is greater than the limiting photon 
flux that is determined by the threshold of the GRB detector. Thus we can determine the 
maximum redshift z max of the burst through 



47rD 2 (z max )(l + z max ) c 



— -Fph.lim, (12) 



where D(z) is the co moving distance of the source to the observer and -F ph lim is the triggering 
threshold of the detector. 

The observed burst rate with isotropic-equivalent energy between E and E + dE is given 

dt ohs dE Jo 1 + z dz v ' 

where -Rgrb is the GRB rate per unit comoving volume per unit comoving time, and V(z) 
is the comoving volume. In a flat cosmology dV(z)/dz = 4tiD 2 (z) dD(z)/dz. 

Nakar et al. (2004) used the flux table of the BATSE 4B Catalog and Band spectrum 
to estimate the distribution of T. They found a lognormal distribution, 



TP(T) = = exp 



where // = 2.15 and a lnT = 0.87 (T = s 



(14) 



Next, we consider the selection effect in redshift identification. Bloom (2003) found that 
there are strong observational biases in ground-based redshift discovery, and suggested to 
use a probability function related to the luminosity distance to minimize this effect. Pl{z) 
is constructed: 

fl, z < z h 

w-xm 1 - *>«■ (15) 

Via K-S test, Bloom constrained z\ < 1.25 and — 3 < L < —1. We here adopt the value of 
Z\ — 1 and allow the parameter L to vary between —3 and —1 case-by-case in order to get a 
better fit. Using the distribution of T and taking into account the redshift selection effect, 
we generalize the distribution of isotropic-equivalent jet energy as 
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Furthermore, since in the ID analysis the integration over redshift hides certain important 
information about the distribution as a function of z, we extend the ID distribution N(E) 
to the 2D distribution N(E, z), 

**->-£^-wr M H ,w ' , ^£ <i7) 



where T max (E,z) is determined by equation (12). 



Finally, we consider the GRB event rate. It has been suggested that GRBs follow 
the star formation rate, because GRBs are probably produced in the final gravitational 
collapse of massive stars (e.g. Woosley 1993; Paczynski 1998). Here we assume that the 
rate of GRBs traces the global star formation history of the universe, Rgrb(z) oc Rsf(z) oc 
Rss(z). Following Porciani & Madau (2001), we employed three different kinds of global 
star formation rate (SFR) model. These three SFR models are similar at z < 1. The 
main difference is at high redshifts: in model 1, the SFR decreases at z > 1.5; in model 
2, the SFR contains the dust extinction effect and remains constant at z > 2; and in 
model 3, the SFR increase at high redshifts after a correction due to a large amount of 
dust extinction. These SFR models are all in an Einstein-de Sitter universe with parameter 
Vt M = 1.0, Vt K = 0.0, Q\ = 0.0, and H = 65 km s~ x Mpc" 1 . We can obtain the SFR in any 
cosmology with parameters M , Q\, fl K and H through the formulation in the appendix 
of Porciani k Madau (2001). 

We assume that all stars with masses M > 3OM explode as Type Ib/Ic SNe, and 
adopt the initial mass function suggested by Madau & Pozzetti (2000). The SN Ib/Ic rate 
can then be estimated by i? S N ib/c = /sn ib/c-RsF = 1-8 x 10 -3 Mq 1 .Rsf- Lamb (1999) showed 
that the observed ratio of the rate of GRBs to the rate of Type Ib/Ic SNe in the observable 
universe is -Rgrb / -^Typc ib/c ~ 10~ 5 . Therefore the GRB event rate can be approximated by 
-Rgrb = /grb-Rsf = lO" 8 M 1 i? SF . 



3. Results 

In order to test the configuration of jets, we firstly perform a theoretical analysis in two 
limits: the rate of GRBs at low and high values of isotropic-equivalent energy. The analysis 
can also predict and test our numerical results. Secondly, we present detailed results of 
numerical calculations, which are well consistent with the theoretical one. 

We used the same parameters as in Perna et al. (2003) and Nakar et al.(2004). The 
triggering threshold is F ph ,iim = 0.424 photons cm" 2 s -1 , and the cosmological parameters 
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are: VL M = 0.3, fi A = 0.7, and H = 71 km s _1 Mpc" 1 



3.1. Theoretical analysis 

The comoving distance in the currently adopted cosmology can be approximated by 
(Wu, Dai, & Liang 2004) 

D(z) ~ . (18) 

V ; H l + 0.29z y ' 

Combined with equation (12), the photon flux in both limits reads 

0.55 photons- crn-.s"^^)" x ^ (19) 

We take the conventional notation Q = Q x ■ 10 X . For a given triggering instrument, Fp^um 
is known. Therefore in the limit of z max < 1 the maximum redshift that a GRB can be 
detected is 

^iwT^-iB' 1 ' 2 - <2o) 

while in the limit of z max > 4, it follows 

*- = 011 fe)" 1 " 53 ^)" 1 ' (21) 

Here we scale the detection threshold to that of the BATSE triggering channel. Assuming 
an empirical T ~ 8 s, the maximum redshift is z oc E 1 ' 2 for low energies E < 10 53 ergs, and 
z oc E for extremely large energies E > 10 54 ergs. In the following analysis, we take the star 
formation rate model 2 as an example. 

In the case of z max < 1, i.e., for GRBs with low isotropic energy of E < 10 53 ergs, the inte- 
gration over redshift z in equation (13) gives 43/o 6 5/GRB,-8(i ? ph,iim/F p B h ^m E )" 3/2 ^53 /2 (^/8 s)~ 3/2 yr" 1 , 
where h e5 = H /Q5 km • s" 1 • Mpc -1 . The differential rate of bursts for the power-law struc- 
tured jet model is 

d2N ( E ) n ( k E \ "V °c \ 2 ^/k E t I F phMl \ ( T \ - 3/2 



cit obs cilog 10 E 



(22) 

Similarly the differential rate for the Gaussian structured jet model is 

d 2 N(E) ( 0c V f ( fph,Um V 3/ V T V 3/2 „ 3 /2 -, ,,,, 



- 9- 



For the uniform jet model, we just present the rate in the case that the distribution of 
isotropic-equivalent energy is an exponential function, 



d 2 N(E) 1/2, o.OU_ 53 ^,5iN * ( -^Ph.lim A ( T \ ' p2 -A_5 3 £ B 3 - 

(24) 

where A = 10~ 53 A_53 erg -1 and £7,51 is the jet energy in units of 10 51 ergs. 

In the case of z max > 4, i.e., E > 10 54 ergs, the integration over redshift z in equation 
(13) gives ~ 1.86 x 10 3 /i^ 5 / G rb -sfi^max) yr _1 . The value of g{z m&x ) ranges from 0.03 to 
0.28. To obtain an estimation, we adopt g(z milx ) ~ 0.1. Then the differential GRB rate for 
the power-law structured jet model is 

d 2 N(E) , „ L . (k E \~ l ( » c \ 2 „ 2/lls , ,j(z 



* 17/4 (55) " (m) E %'J°™--* gJ ^r E * ^- < 25 > 



dt ohs dlog w E 

Similarly, the rate for the Gaussian structured jet model is 

d 2 N(E) ^ noK , 4 , 9(z m&x )f9, 



dt ohs dlog 10 E 



-0.85^ GRB ,_ 8 ^(^ 2 yr-. (26) 



At last, for the uniform jet model, if the distribution of jet isotropic-equivalent energy con- 
forms to an exponential distribution, then the differential burst rate is 

~7T F - 60.4Z&A ,: i ('"" IA - A ' •)/^.A;h>, >%^4/ J < A ^ yr" 1 . (27) 

dt ohs dlog 10 E 0.1 

Figure 1 shows differential rates of GRBs determined by equation (13). Different lines 
correspond to the predicted distributions derived from different possible models. It can 
be seen that the above analytical expressions describe the rate quite well at low and high 
limits of E. The predicted GRB rates for these models differ from each other significantly, 
especially at extremely low and high energy limits. Such evident differences enable us to find 
the best model for the jet structure, if the sample contains enough bursts. Therefore the ID 
comparison between the predicted distribution of E and the observed sample is available, 
although the integration over redshift conceals much information. 



3.2. Numerical analysis 



To compare the theoretical distributions derived in §2, we consider a sample of GRBs 
whose fluences and redshifts are measured. Such a sample including 41 bursts are listed in 
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Table 1. This table does not contain GRBs such as GRB 980329, 980519, 000630, 020331, 
030115 and 040511, because of the large uncertainties of their redshifts. We rule out 2 GRBs 
in this table: GRB 980425 and 031203, because they form a peculiar subclass characterized 
by their unusually low isotropic gamma-ray energy releases and other unusual properties. 
For example, these two GRBs violate the P 7! i S o — Ppeak relation and the luminosity-spectral 
lag relation (Sazonov, Lutovinov, & Sunyaev 2004). Furthermore, Yamazaki, Yonetoku and 
Nakamura (2003) found that the observed unusual properties of the prompt emission of GRB 
980425 could be explained by using an off-axis jet model of GRBs, such as the extremely low 
isotropic equivalent 7-ray energy, the low peak energy, and the high fluence ratio, and the long 
spectral lag. Ramirez- Ruiz et al. (2004) also argued that the observed data for GRB 031203 
are more consistent with a GRB seen at an angle of about twice the opening angle of the 
central jet. Since we do not consider the off-axis case in our model analysis, we rule out these 
two GRBs in the numerical analysis. Therefore, we finally have 39 bursts available in our 
analysis. We have calibrated the isotropic-equivalent energy release in a fixed cosmological 
bandpass, 0.1 keV-10 MeV, using the cosmological k-correction (Bloom, Frail, & Sari 2001). 
The calorimetric isotropic-equivalent energy is given by E = 47tDl5'[o.i-io 4 ]/(1 + z )- 

We calculate the ID distribution N(E) under three jet structure scenarios, and perform 
K-S test to assess the compatibility of the theoretical distributions with the observed data. 
We describe the significance level for the result of K-S test as Pk- S - Small values of Pk_ s show 
that the theoretical distributions and the observed data are significantly different (Press et 
al. 1997). In the analysis below, we take the star formation rate model 2 and L = —1.0 as 
an example. We also perform the test for other star formation rate models as well as other 
values of L, and ultimately obtain similar results. 

Figure 2 shows that the predicted distribution peaks approximately at the same energy 
for different SFR models. Three SFR models result in slight differences in the predicted 
distribution if we adopt other parameters with the same values. In Figure 3 the 2D grey 
contours show the confidence level (for the null hypothesis that the observed data are drawn 
from this model) as a function of fog and 9 C . We find that the predicted distribution is 
consistent with the observations at the confidence level of > 40% while the parameters fc^ 
and 6 C are changed in a wide range. For other star formation rate models, the power-law 
structured jet model also agrees with the observational data reasonably well. 

However, the distribution derived from the Gaussian jet model does not agree with 
the observations even if the parameter 6 C is changed in a wide range. Figure 4 shows the 
theoretical distribution of GRB isotropic energy for different star formation rate models. The 
theoretical distribution peaks at higher energy compared with the observed sample. Figure 5 
shows the predicted and observed cumulative distributions of GRB isotropic energy for three 
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different values of L, assuming the star formation rate model 2. Different star formation 
rate models and different values of L only lead to a slight difference in the distribution. 
The Gaussian jet model predicts that most GRBs are detected within its central core where 
the isotropic energy is large and nearly constant. This is because the energy decreases 
exponentially at the wing and the GRB cannot be triggered by the detectors. However, 
we note that if we introduce some dispersion on the value of the central angle 9 C , this 
inconsistency may be reduced (Zhang et al. 2004; Dai & Zhang 2005). 

In the uniform jet model, we have considered three cases in which f(E) is a power- 
law, exponential, or Gaussian function. By comparing the theoretical calculation and the 
observational data using the K-S test, we obtain a confidence level for the hypothesis that 
the observational data is drawn from the theoretical distribution. In the case of a power- 
law distribution, the result of the K-S test depends mainly on the value of the index (. 
Our numerical results show that when ( varies between -0.99 and -0.79, the theoretical 
distribution is compatible with the observed data at the confidence level of > 40%. This 
confidence reaches a maximum of 90% for ( = —0.91. We show in Figure 6 the predicted 
and observed distributions of GRB isotropic energy for three different values of (. Figure 
7 presents the predicted and observed cumulative distributions of GRB isotropic energy for 
three different values of (. 

In the case of exponential and Gaussian distributions, the results of the K-S test also 
depend on the parameters. But the theoretical distributions are consistent with the obser- 
vations at the confidence level of < 5%, even though the parameters are changed greatly. 
Figures 8 and 9 show the differential and cumulative distributions of GRB isotropic energy 
respectively. Despite the variation of the parameter A, the figures show a poor consistency 
between the theoretical prediction and the observed data. Figures 10 and 11 exhibit similar 
results for the case that f(E) follows a Gaussian function. 

We next consider the 2D distribution N(E, z). The distribution of GRBs is a function of 
both isotropic equivalent energy and redshift. In the ID analysis, the information contained 
in the redshift space is concealed by the integration over redshift. In order to explore the 
overall information about the distribution of GRBs, we extend the distribution from one 
dimension to two dimensions, and compare it with the observed sample (see Figs. 12-16). 
We performed 2D K-S test to check the consistency of the predicted 2D distribution with 
the observations. For the sample of 39 bursts, a 2D K-S test shows that for the power-law 
structured jet model, the predicted distribution is consistent with the observations at the 
confidence level of > 40% while the parameters k E and 6 C are changed in a wide range. 
When adopting the value of parameters k E = 2.2 and 6 C = 0.02, we obtain the best fit, at 
the confidence level of 58%. For the uniform jet model in the power-law case, the model 



- 12 - 



with SFR2, L = —1, and ( = —0.92 corresponds to a best fit, at the confidence level of 
57%. When ( varies between -0.97 and -0.85, the predicted distribution is compatible with 
the observed data at the confidence level of > 40%. In the exponential and Gaussian case, 
we found that the hypothesis that the data are drawn from these three models is completely 
rejected. This conclusion is unchanged for different values of the parameters. 

However, this conclusion is tentative, because the sample suffers from various observa- 
tional biases, such as GRB detection, localization, and especially the selection effect in the 
identification of redshift (Bloom 2003). Nakar et al. (2004) suggested that the selection 
effect in redshift can be minimized by testing the 9 distribution for a given redshift. But the 
size of the sample is greatly reduced in this method. With the cumulation of the sample in 
the future, this method might be used to test the configuration of jets in GRBs. 

4. Discussion and Conclusion 

In this paper, we have worked out the theoretical distribution of burst isotropic-equivalent 
energy in ID N(E) = dN/dE and in 2D N(E,z) = dN 2 / ' dEdz for several jet structure 
models. The figures and analyses above show that theoretically the GRB rates predicted 
by different models differ from each other greatly. It is possible for us to find the true 
model by the method shown here if the sample is large enough. Based on the theoretical 
analysis in §3.2, we carried out numerical calculations and compared N(E) = dN/dE and 
N(E, z) = dN 2 /dEdz with the the observed sample of GRBs with known fluences and red- 
shifts. Via K-S test, we found that the power-law structured jet model and the uniform jet 
model can be consistent with the currently observed data at a confidence level of ^ 40% 
under certain circumstances. 

One advantage of our work over Perna et al. (2003) and Nakar et al. (2004) lies in the 
fact that the sample of bursts with measured redshifts and fluences is larger than the one 
with known redshifts and modelled angles. The sample of bursts with known redshift and 
modelled angle is too small (currently 16 GRBs). According to the relation tj oc 9^ 3 (Sari 
et al. 1999), if the jet opening angle 9j is small, the jet break time tj probably is so early 
that there is only the upper limit on it. In the contrary, if 9j is large, the late-time optical 
transient may be too dim to be detected. As a result, this selection effect in detecting tj 
leads to a small sample of bursts with measured jet break times. However, the sample of 
bursts with measured fluences and redshifts does not suffer from such selection effect. This 
larger sample helps us to draw a stronger conclusion. 

However, it is worth noting that all results above are obtained under the following 
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assumptions. First, the isotropic-equivalent energy release in gamma-ray band is assumed to 
be proportional to the total isotropic-equivalent energy of the bursts. Second, we consider 
the simplifying assumption that the theoretical model is strictly consistent with the actual 
jet structure. These theoretical models predict a smooth, broken power-law light curve. 
However, recent observations revealed variabilities in some afterglow light curves such as in 
GRB 021004 (Fox et al. 2003). One leading model named "patchy shell" model attributes 
this fluctuation in light curves to random angular fluctuation of the energy per solid angle 
(Kumar & Piran 2000; Nakar, Piran, & Granot 2003). In this paper, we consider a regular 
energy distribution in the jet rather than with some random angular fluctuations. This 
simplifying assumption may result in an underestimate of the number of energetic events. 

Third, the current sample of 39 bursts with measured redshifts suffers from numerous 
selection effects, especially in optical afterglow detections and redshift identifications. To 
reduce these effects, we used a strategy suggested by Bloom (2003). On the other hand this 
method might introduce a systematic error into the result. To reach a robust conclusion, 
Nakar et al. (2004) put forward an approach to minimize the selection effect in z, i.e., to 
test the distribution for a given z. In our work, to avoid this selection effect and reduce the 
dependence on SFR model, we need to derive the distribution of isotropic-equivalent energy 
for a fixed redshift N(E)\ Z . However, this method requires a large sample of bursts with 
measured fluences and with the same redshift, while there are only a few bursts are detected 
with the same redshift. With a larger sample this method shown here can be expected to give 
a stronger conclusion. With the advent of the Swift era, the sample of GRBs with measured 
fluence and redshift is anticipated to be much larger than the current one. The large sample 
of bursts detected by the same detector will be having the same parameters such as the 
detector sensitivity, thus reducing the uncertainty and reaching a stronger conclusion. 

We thank E. W. Liang, D. Xu and Z. P. Jin for their helpful discussions. This work 
was supported by the National Natural Science Foundation of China (grants 10233010 and 
10221001) and the Ministry of Science and Technology of China (NKBRSF G19990754). 
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Table 1. Compilation of Spectra and Energetics Input Data 



GRB a 


z 




Bandpass 


K/3] c 


E° hs d References 






[1CT 6 org cm- 2 ] 


[keV] 




[kev] (z, S 7 = S, tjet = t, n, a, (3, E v ) 



970228 


0.6950 


11.00 


970508 


0.8349 


1.80 


970828 


0.9578 


96.00 


971214 


3.4180 


8.80 


980326 


1.0000 


0.75 


980425 


0.0085 


3.87 


980613 


1.0969 


1.00 


980703 


0.9662 


22.60 


990123 


1.6004 


300.00 


990506 


1.3066 


194.00 


990510 


1.6187 


19.00 


990705 


0.8424 


75.00 


990712 


0.4331 


11.00 


991208 


0.7055 


100.00 


991216 


1.0200 


194.00 


000131 


4.5000 


35.10 


000210 


0.8463 


61.00 


000214 


0.4200 


1.42 


000301c 


2.0335 


2.00 


000418 


1.1182 


20.00 


000911 


1.0585 


230.00 


000926 


2.0369 


6.20 


010222 


1.4769 


120.00 


010921 


0.4509 


18.42 


011121 


0.3620 


24.00 


011211 


2.1400 


5.00 


020124 


3.1980 


8.10 


020405 


0.6899 


74.00 


020813 


1.2550 


102.00 


021004 


2.3351 


2.55 


021211 


1.0060 


3.53 


030226 


1.9860 


5.61 


030323 


3.3718 


1.23 


030328 


1.5200 


36.95 


030329 


0.1685 


110.00 


031203 


0.1055 


1.20 



40, 700 
40, 700 

20, 2000 
40, 700 
40, 700 

20, 2000 
40, 700 

20, 2000 
40, 700 

20, 2000 
40, 700 
40, 700 
2, 700 

25, 1000 
20, 2000 

26, 1800 
40, 700 
40, 700 

150, 1000 
15, 1000 
15, 8000 
25, 100 
2, 700 
2, 400 
25, 100 
40, 700 
2, 400 
15, 2000 
30, 400 
2, 400 
2, 400 
2, 400 
2, 400 
2, 400 
30, 400 
20, 2000 



-1.54, 


-2.50 


115 


z: 


1, S: 2, a: 2, (3: 


2, E p : 


2 




-1.71, 


-2.20 


79 


z: 


3, 


S: 2, a: 2, (3: 


2, E p : 


2 




-0.70, 


-2.07 


298 


z: 


4, 


S: 5, a: 6, /3: 


5, E p : 


5 




-0.76, 


-2.70 


155 


z: 


7, 


9: 


2, a: 2, f3: 


2, E p : 


2 




-1.23, 


-2.48 


47 


z: 


8, S: 


2, a: 2, (3: 


2, E p : 


9 




-1.27, - 


2.30 * 


118 


z: 


10, 


s 


5, a: 5, E 


P • 5 






-1.43, 


-2.70 


93 


z: 


11, 


s 


2, a: 2, f3 


2, E p 


2 




-1.31, 


-2.40 


254 


z: 


12, 


s 


5, a: 6, (3 


5, E p 


5 




-0.89, 


-2.45 


781 


z: 


13, 


s 


2, a: 2, f3 


2, E p 


2 




-1.37, 


-2.15 


283 


z: 


14, 


s 


5, a: 6, (3 


5, E p 


5 




-1.23, 


-2.70 


163 


z: 


15, 


s 


2, a: 2, f3 


2, Ep 


2 




-1.05, 


-2.20 


189 


z: 


16, 


s 


2, a: 2, (3 


2, Ep 


2 




-1.88, 


-2.48 


65 


z: 


15, 


s 


17, a: 2, (3: 2, E p : 2 










z: 


18, 


s 


19 








-1.23, 


-2.18 


318 


z: 


20, 


s 


5, a: 6, f3: 5, E p 


5 




-1.20, 


-2.40 


163 


z: 


21, 


s 


21, a: 21, /3: 21, 


Ep\ 


21 








z: 


22, 


s 


22 








-1.62, 


-2.10 


> 82 


z: 


23, 


s 


2, a: 2, f3: 2, E p 


2 










z: 


24, 


s 


25, 














z: 


14, 


s 


26 








-1.11, 


-2.32 


579 


z: 


27, 


s 


27, a: 27, (3: 27, E p : 


27 








z: 


28, 


s 


29 








-1.35, 


-1.64 


> 358 


z: 


30, 


s 


31, a: 2, (3: 2, E p : 2 




-1.55, 


-2.30 


89 


z: 


32, 


s 


33, a: 33, f3: 34, 


Ep\ 


33 


-1.42, - 


2.30 * 


217 


z: 


35, 


s 


36, a: 6, E p : 6 






-0.84, - 


2.30 * 


59 


z: 


37, 


s 


37, a: 6, E v : 6 






-0.79, 


-2.30 


120 


z: 


38, 


s 


33, a: 33, (3: 34, E p : 


33 


0.00, 


-1.87 


192 


z: 


39, 


s 


39, a: 39, f3: 39, E p : 


40 


-1.05, 


-2.30 


212 


z: 


41, 


s 


42, a: 42, f3: 42, E p : 


42 


-1.01, - 


2.30 * 


80 


z: 


43, 


s 


33, a: 33, E p : 33 




-0.80, 


-2.37 


47 


z: 


44, 


s 


45, a: 45, f3: 45, E p : 


45 


-0.89, 


-2.30 


97 


z: 


46, 


s 


45, a: 45, f3: 34, E p : 


45 


-1.62, - 


2.30 * 




z: 


47, 


s 


33, a: 33 








-1.14, 


-2.09 


126 


z: 


48, 


s 


33, a: 33, f3: 33, E p : 


33 


-1.26, 


-2.28 


68 


z: 


49, 


s 


50, a: 50, (3: 50, 


Ep\ 


50 


1.00 *, 


-2.30 * 


> 190 


z: 


51, 


s 


52, E p : 53 







Table 1 — Continued 



GRB a 



z 



[10 ^ erg cm 2 




E° hs d References 

[kev] (z, 5 7 = S, tjet = t, n, a, f3, E p ) 



040924 

041006 

050408 

050525a 

050603 



0.8590 
0.7160 
1.2357 
0.6060 
2.2810 



2.73 
19.90 

1.90 
20.00 
34.10 



20, 500 
25, 100 
30, 400 
15, 350 
20, 3000 



1.17, -2.30 * 
1.37, -2.30 * 
1.98, -2.30 * 
1.00, -2.30 * 
-0.79, -2.15 



67 
63 
20 
79 
349 



z: 54, S: 55, a: 56, E p : 55 

z: 57, S: 58, a: 56, E p : 56 

z: 59, S: 59, a: 60, E p : 60 

z: 61, S: 61, a: 61, E p : 61 



z: 62, S: 63, a: 63, /3: 63, E p : 63 



a Upper/lower limit data are indicated with < and > respectively. References arc given in order for redshift ("z"), 
fiuence ("S"), low energy band spectral slope ("a"), high energy band spectral slope ("/?"), and spectral peak energy 



b GRB fiuence 5 7 calculated in the observed bandpass [ei, e2] keV. 

c Low energy "Band" spectral slope a and high energy "Band" spectral slope (3. When j3 is reported in the literature 
but a is not, wc set a = —1.00 (marked with *). Following Attcia (2003), when a is reported in the literature and (3 
is not, we fix /3=-2.30 (marked with *). When both a and (3 are not reported in the literature we assume o=-1.00 and 
/3=-2.30. 

d Observed spectral peak energy E* s = E° hs (2 + a). When E° hs is not reported we fix B* s = 250 keV. When E° hs 
is reported only with a limit, we adopt E° hs as the value of the limit. 
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Fig. 1. — The observed differential rate of GRBs as a function of E for 5 kinds of jet model. 
The parameter T is chosen to be 8 s for all GRBs. Other parameters are assumed: L = —1, 
Ej = 10 51 ergs, /grb = 10~ 8 Mq 1 , and -F phi i im = 0.424 photons • cm" 2 ■ s -1 . For power- law 
structured jet models, = 2.0 and 9 C = 0.01. For the Gaussian structured jet model, 
8 C = 0.05. For the uniform jet model A = 0.7 x 10~ 53 ergs -1 (the exponential case); ( = —1.1 
(the power-law case); E = 0.7 x 10 53 ergs, a = 1.4 x 10 53 ergs (the Gaussian case). Here we 
take the star formation rate model 2 as an example. 
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Fig. 2. — Theoretical distribution of the isotropic-equivalent energy in the power-law struc- 
tured jet model (lines), compared with the observed distribution from a sample of 39 GRBs 
with known redshifts detected so far (histogram). Different lines correspond to different star 
formation rate models. The parameters: &e = 2.3 and 9 C = 0.02. Here we take L = — 1 as 
an example. 
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Fig. 3. — The Gray scale map of the K-S test result for different values of Ue and 6 C . Every 
point in this plot corresponds to a set of parameters for the PLSJ model and the color of 
the point represents how well the PLSJ model with such parameters is compared with the 
observed data. Here we take the star formation rate model 2 and L = — 1 as an example. 
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Fig. 4. — Theoretical distribution of the isotropic-equivalent energy in the Gaussian struc- 
tured jet model (lines), compared with the observed distribution from a sample of 39 GRBs 
with known ffuences and redshifts detected so far (histogram). Different lines correspond to 
different star formation rate models. The parameter 6 C = 0.01. Here we take L = — 1 as an 
example. 
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Fig. 5. — Cumulative distribution for the Gaussian structured jet model. Three values of 
L, the index in the redshift-identification selection effect, are considered. The cumulative 
histogram is plotted from the sample. The parameter 6 C = 0.01. Here we take star formation 
rate model 2 as an example. 
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Fig. 6. — Theoretical distribution of the isotropic-equivalent energy in the uniform jet model 
in the case that the distribution of isotropic-equivalent energy is a power-law function (lines), 
compared with the observed distribution from a sample of 39 GRBs with known fluences 
and redshifts detected so far (histogram). Different lines correspond to different values of 
the power-law index (. Here we take the star formation rate model 2 and L = — 1 as an 
example. 
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Fig. 7. — Cumulative distribution for the uniform jet model in the case that the distribution 
of isotropic-equivalent energy is a power-law function. Different lines correspond to the 
different values of the power-law index £. The cumulative histogram is plotted from the 
sample. Here we take the star formation rate model 2 and L = — 1 as an example. 
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Fig. 8. — Theoretical distribution of the isotropic-equivalent energy in the uniform jet model 
in the case that the distribution of isotropic-equivalent energy is an exponential function 
(lines), compared with the observed distribution from a sample of 39 GRBs with known 
fluences and redshifts detected so far (histogram). Different lines correspond to different 
values of the exponential parameter A. Here we take the star formation rate model 2 and 
L = — 1 as an example. 
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Fig. 9. — Cumulative distribution for the uniform jet model in the case that the distribution 
of isotropic-equivalent energy is an exponential function. Three values of A are considered. 
The cumulative histogram is plotted from the sample. Here we take the star formation rate 
model 2 and L = — 1 as an example. 
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Fig. 10. — Theoretical distribution of the isotropic-equivalent energy in the uniform jet model 
in the case that the distribution of isotropic-equivalent energy is a Gaussian function (lines), 
compared with the observed distribution from a sample of 39 GRBs with known fluences 
and redshifts detected so far (histogram). Different lines correspond to different values of 
mean energy E (in the units of ergs). We adopt the scatter in the Gaussian function as 
a = 1.7 x 10 53 ergs. Here we take the star formation rate model 2 and L = — 1 as an 
example. 
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Fig. 11. — Cumulative distribution for the uniform jet model in the case that the distribution 
of isotropic-equivalent energy is a Gaussian function. Different lines correspond to different 
values of mean energy E (in the units of ergs). We adopt the scatter in the Gaussian function 
as cr = 1.7 x 10 53 ergs. The cumulative histogram is plotted from the sample. Here we take 
the star formation rate model 2 and L = — 1 as an example. 
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Fig. 12. — The 2D distribution density, N(E, z), for the power-law structured jet model and 
star formation rate model 2. The parameters are assumed: L = —1, )ze = 2.2, and 8 C = 0.02. 
The circles denote 39 bursts with known fluences and redshifts detected so far. 
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Fig. 13. — The 2D distribution density, N(E, z), for the Gaussian structured jet model, star 
formation rate model 2. The parameters are assumed: L = —1 and 9 C = 0.05. The circles 
denote 39 bursts with known ffuences and redshifts detected so far. 
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Fig. 14. — The 2D distribution density, N(E,z), for the uniform jet model (the power-law 
case) and star formation rate model 2. The parameters are assumed: L = —1 and £ = —0.9. 
The circles denote 39 bursts with known fluences and redshifts detected so far. 
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Fig. 15. — The 2D distribution density, N(E,z), for the uniform jet model (the exponential 
case) and star formation rate model 2. The parameters are assumed: L — — 1 and A = 
0.7 x 1CT 53 ergs -1 . The circles denote 39 bursts with known fluences and redshifts detected 
so far. 
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Fig. 16. — The 2D distribution density, N(E,z), for the uniform jet model (the Gaussian 
case) and star formation rate model 2. The parameters are assumed: L — — 1 and E = 
0.7 x 10 53 ergs, and o = 1.0 x 10 53 ergs. The circles denote 39 bursts with known ffuences 
and redshifts detected so far. 



